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Gamow- Teller (GT) strength distributions of Mg isotopes are investigated within a framework of 
the deformed quasi-particle random phase approximation(DQRPA). We found that the N=20 shell 
closure in ^*~^^Mg was broken by the prolate shape deformation originating from the /p-intruder 
states. The shell closure breaking gives rise to a shift of low-lying GT excited states into high-lying 
states. Discussions regarding the shell evolution trend of single particle states around N=20 nuclei 
are also presented with the comparison to other approaches. 



PACS numbers: 
Keywords: 



23.40. He, 21.60.Jz, 26.50.-|-x 



The structure of nuclei along the valley of stability is 
well understood in terms of the conventional magic num- 
bers in the nuclear shell model. In a region far from the 
stability, experimental studies have shown the breaking 
of the N=20 magic number in N=20 nuclei, such as "^^Na, 
■^^Mg, and '^■^Al nuclei located in so called the island of 
inversion [1 - 3] . The breaking comes from the burrowing 
of /7/2 and P3/2 states below ^3/2 state f3|, which may 
originate from the monopole component of the nucleon- 
nucleon interaction by the tensor force Q ■ The breaking 
of magic numbers in the inversion island is thought as 
one of the interesting properties of the exotic nuclei near 
drip line. 

The neutron-rich nuclei around N=20 are thought to 
be strongly deformed because of the reduction of the 
N=20 shell gap. The deformation in nuclei becomes more 
important than any other periods with the recent devel- 
opment of rare isotope accelerator facilities, from which 
one may perform lots of challenging experiments related 
to the exotic nuclei. Although they decay fractions of a 
second, their existence is imprinted on the nuclear abun- 
dances of stars through the successive nuclear reactions 
in the cosmos at the explosive stage of stellar evolution, 
such as supernovae (SNe) explosion 0. 

For example, the nuclear reactions induced by the neu- 
trinos emitted from the SNe are treated as important 
input data for the neutrino(z^)-process Since the 

neutrinos emitted from a proto-neutron star may have 
tens of MeV energy high enough to excite the deformed 
nuclei, one needs to understand more precisely the high- 
lying (HL) excited states beyond one nucleon threshold. 
Among them, the GT states are of great importance be- 
cause most of the charge exchange reactions (CEXRs) arc 
dominated by the GT transition. These HL GT excited 
states are in close association with the nuclear structure. 

One of typical instances is the GT quenching problem. 
It says that the difference of total running sums S± for 



GT (±) transitions, (5_ —5+), is usually quenched com- 
pared to the Ikeda sum rule (ISR), {S- -S+) = 3{N-Z). 
But, recent experimental data on HL GT states deduced 
by more energetic projectiles shed a new light on the GT 
states located above one nucleon threshold, whose con- 
tributions enable us to explain the quenching problem 
through the multi-particle and multi-hole configuration 
mixing 

Conventional approach to understand the nuclear 
structure is based on the spherical symmetry. But, in 
order to describe neutron-rich nuclei and their relevant 
nuclear reactions occurred in the nuclear processes, one 
needs to develop theoretical frameworks including explic- 
itly the deformation [ll|, [T^. Ref. 11 1 used an effec- 
tive separable force to the deformed quasi-particle ran- 
dom phase approximation (DQRPA). Realistic two-body 
interaction was firstly exploited at Ref. [l3], but only 
with neutron-neutron (nn) and proton-proton (pp) pair- 
ing correlations which have only isospin T = 1 and J = 
interaction. But, to properly describe the deformed nu- 
clei, the T= and J = 1 pairing should be also taken 
into account because the J = 1 pairing is believed as one 
of the reasons leading to the nuclear deformation. 

In this paper, we apply our previous DQRPA [H, [l3l 
to the nuclei in the N=20 inversion island. Starting from 
the deformed Wood Saxon potential [l^ , we transform a 
physical state given by the diagonalization of total Hamil- 
tonian in the Nilsson basis into the spherical basis, in 
which one can perform more easily theoretical calcula- 
tions. Then deformed BCS equation can be solved by 
using ab initio Brueckner G-matrix represented in the 
deformed basis for the nucleon-nucleon interaction inside 
nuclei as follows. 

With the expansion coefficient Ba from the deformed 
to the spherical basis, lafla >= J2a^a W^a >, the 
pairing potential Ap is calculated as 
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FIG. 1: (Color online) Neutron single particle energies of 
N=20 isotones as a function of the proton number with the 
deformation parameter /32 from the relativistic mean field 
(RMF) (upper panel) and from B(E2) [13, [Hi (lower panel) . 



troduced for the transformation to the deformed basis of 
G-matrix. Here K' , which is a projection number of the 
total angular momentum J onto the z axis, is selected 
isT' = at the BCS stage because we consider the pair- 
ings of the quasi-particles at a and a states. In order 
to renormalize the G-matrix, strength parameters, g^air^ 
9pair ^^'^ ffpair ^-^^ multiplied to the G-matrix [Tsj . 

The transition amplitudes from the ground state of 
an initial nucleus to the excited state are expressed by 

< l{K),m\j3^\ QRPA > 

= ^ I^aa'-p^PP'-pf, < Cta"ppa\(7K\fiP"npp > 

x[Upaa"Vnl3l3"Xaa"l3l3".K + Vpaa"Unl3l3"Yaa" 13/3" mJ."^) 

where Af is the nomalization factor and | QRPA > de- 
notes the correlated QRPA ground state in the intrinsic 
frame. 



FIG. 2: (Color online) Neutron single particle energies of Mg 
isotopes as a function of the neutron number with the defor- 
mation parameter /32 from the relativistic mean field (RMF) 
(upper panel). In lower panel, the deformation parameter 
l^^mp j^j. 26j^g ^32j^g and ='''Mg are from the B(E2) data in 
Refs. [TH and [3, respectively. 

To compare to the experimental data, the GT(±) 
strength functions 

^gtM = E I < UK),m\0^\\ QRPA > p (3) 

and their running sums Sqj, = Y,raBQj,{rn) are eval- 
uated. The single particle states are used up to AItuj 
in the spherical limit. Since the GT strength distribu- 
tions turn out to be sensitive on the deformation pa- 
rameter, [IBl, we took 02 values from both the rel- 
ativistic mean field (RMF) [16] and the empirical data, 
/32 = (47r/3Zi?2)[S(^2 ^)/e'Y/^{RQ = l.2A^'^)^^. 

In order to compare our single particle state energies 
(SPSEs) of N = 20 nuclei to those by other nuclear mod- 
els, which show the shell evolution, we plotted SPSEs in 
Figs. 1 and 2, whose deformation parameters are tabu- 
lated in Table I. 
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TABLE I: Deformati on p arameters /32 of N=20 isotones from 
RMF and B(E2)[ia[ii| for N=20 isotones. 



A 



jRMF /3=-P(B(E2)) 



30Ne 
3*Si 

36s 

3«Ar 

Ca 



40 



0. 
0. 
0. 
0. 
0. 
0. 



0.580 
0.473 
0.179 
0.168 
0.163 
0.123 



TABLE II: Deformation parameters and /JJ'"", and em- 

pirical(theoretical) pairing gaps Af;^(A^^") for Mg isotopes. 
Tiie particle-particle(particle-hole) strength parameters are 
taken as gpp = 1.0{gph ~ 1.0). 



A 


qRMF 
P2 


nemp 
P2 


AS™ (AfJ 




30 Mg 


0.17 


0.431 


2.094(2.095) 


1.700(1.709) 


^^Mg 





0.473 


2.391(2.42) 


1.645(1.649) 


^*Mg 


0.192 


0.54 


3.29(3.308) 


1.895(1.912) 



The level density of /p shell in the upper panel in Fig.l 
is increased with the decrease of proton number similarly 
to the result of Ref. ^ which exploited HF with SkM 
force (see Fig. 2 at Ref. [1^) in the QRPA scheme to 
study the first 2+ states in ^^Mg and ^"Ne. This trend 
becomes manifest if we use the empirical deformation 
parameter (lower panel) . Of course, the level density may 
lead to the anomaly of the B(E2) values as shown in Ref. 
[Tsl l. In particular, N=20 shell gap is explicitly broken 
with the decrease of proton number by the intruder of 
fi/2 state below state. Recent calculation by the 
shell model (see Fig. 3 (a) at Ref. f2l'|) also shows such 
a trend, although the effect itself was smaller than ours. 

In the following, we discuss the GT strength distribu- 
tions in Mg isotopes, in which ■^^Mg is thought to locate 
in the island inversion. In Table II, /Jg™'' and 

relevant pairing gaps are tabulated. Changes of SPSEs 
by the prolate deformation are presented in Fig. 2. The 
most remarkable point in the lower panel, which used 
/Jj™^ accounting for the large B(E2) strength, is the dis- 
appearance of the magic number N = 20 for N = 16 ^ 
22 nuclei. 

These shell evolution phenomena by the intruder are 
intimately related to the following facts. The SPSEs 
adopted from the deformed Woods Saxon potential nat- 
urally depend on the parameter /32. The deformation 
of nuclei may be conjectured to come from macroscopic 
phenomena, for example, the core polarization, the high 
spin states and so on. Microscopic reasons may be traced 
to the tensor force in nucleon-nucleon interaction, which 
is known to account for the shell evolution according to 
the recent shell model calculations [2l|, . For example. 



T = 0, J = 1 pairing, which is associated with the So 
tensor force, may lead to the deformation compared to 
the T = 1, J = pairing. Therefore, the deformation pa- 
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FIG. 3: (Color online) Gamow- Teller strength distribution 
B(GT~) on Mg isotopes as a function of the excitation en- 
ergy Eex w.r.t. the ground state of '^"Mg ^''^Mg with the 



deformation parameter 



in the upper panel in Fig. 2. 



rameter adopted in this work may include implicitly and 
effectively such effects, because the single particle states 
from the deformed Wood Saxon potential show strong 
dependence on the /32, as shown in Figs. 1 and 2. 

Now we discuss the GT strengths on Mg isotopes whose 
SPSEs heavily depend on the neutron numbers as shown 
in Fig. 2. Since the nuclei considered here are expected 
to have large energy gaps between proton and neutron 
spaces, similarly to the neutron-rich nuclei of importance 
in the r-process, we consider only the nn and pp pairing 
correlation. But for the p-process nuclei, the np pairing 
could be more important than the neutron-rich nuclei 
because of the adjacent energy gaps of protons and neu- 
trons. The calculations for the neutron-deficient nuclei 
in p-process are in progress. 

In Figs. 3 and 4, the GT strength distributions, 
B(GT^), on Mg isotopes are presented as a function 
of the excitation energy E^x w.r.t. the ground state of 
30Mg -34j^g for [Ifi] and P^"^'' values from B(E2) 

data. The GT strength distributions are widely scattered 
into the HL states owing to the deformation. Actually, 
these HL GT states are intimately associated with the 
2p-2h contributions which result from the wide smearing 
of some physical states around the Fermi surface by the 
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FIG. 4: (Color online) Same as Fig. 3 but with the deforma- 

3er. 
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tion parameter /32'"'' in the lower panels in Fig. 2. 



prolate deformation (T^, [3 • 

One more interesting point is that some low-lying GT 
states disappeared if we use empirical /32 values. This 
can be explained if we recollect the disappearance of the 
magic number iV = 20 in Fig. 2. For detailed analysis, 
in Fig. 5, we identify main GT transitions on '^^Mg for 
p^mp _ 4^3^ which are the forward amplitudes X of 
three excitation energies having relatively larger B(GT) 
values in Fig. 4 (b). Main transition in K = comes 
from (2025/2) state, while major transitions in K = 1 are 
(2113/2)-j-(2111/2) and (3213/2)-j>(3301/2). An impor- 
tant point to be noticed in Fig. 5 is that the (3301/2) state 
stemming from the /7/2 state may contribute to the GT 
transition at 10.6 MeV. The (3301/2) state is generated 
below Fermi surface energy by the /j»-intruder states. 

In summary, to describe single particle states in de- 
formed basis, we exploited the deformed axially symmet- 
ric Woods-Saxon potential, and performed the deformed 
BCS and deformed QRPA with the realistic two-body in- 
teraction recalculated in deformed basis from ab initio G- 



matrix based on Bonn potential. Results of the Gamow- 
Teller strength for Mg isotopes show that the deforma- 
tion effect leads to a fragmentation of the GT strength 
into HL states and predicts the HL GT excitations to be 
explored by higher energy projectiles, whose excitations 
may affect the nuclear process in the cosmos through 
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FIG. 5: (Color online) Main transitions in the forward am- 
plitudes X for three dominant GT states w.r.t. the ground 
state of "^^Mg as a function of the configuration state for the 
case of Fig. 4 (b). 



the GT transition. They are shown to result from the 
wide smearing by the increase of the Fermi surface en- 
ergy due to the prolate deformation. This work was sup- 
ported by the National Research Foundation of Korea 
(2012R1A1A3009733, 2012-009733 and 2011-0015467). 
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